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Abstract. We have detected the (0, l lc ,0) J = 3-2 and 
(0,0,0) J = 8-7 transitions of HCN towards the S star 
X Cygni. The excitation requirements of these transitions 
are too high to be satisfied in the outer envelope of the 
star, and the emission must originate within < 20 stellar 
radii, ie. the molecule must be forming close to the star. 
This conclusion is supported by a model for AGB stars in 
which molecules including HCN form in a shocked wind 
close to the stellar surface. 
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1. Introduction 

Recent detections of warm carbon bearing molecules in 
different O-rich AGB stars suggest that these carbon 
species have to form in the deep layers of the stellar wind. 
Photochemical models (Willacy & Millar 1997) of several 
O-rich stars, which rely on injection of certain molecules 
to generate a carbon-rich chemistry at large radii, succeed 
in reproducing the observed values of certain molecules 
but fail to reproduce some molecular abundances, in par- 
ticular that of HCN. HCN is observed in the envelope 
of several O- and S-type stars at abundances of > 10 -6 
(Bieging & Latter 1994 (BL94); Bujarrabal et al. 1994; 
Sopka et al. 1989). 

Duari et al. (1999) (DCW99) considered whether a 
shocked wind model could reproduce HCN abundances for 
O-rich objects. Local thermodynamic equilibrium (LTE) 
models fail to reproduce the observed abundances by sev- 
eral orders of magnitude (BL94) . An investigation of non- 
equilibrium chemistry of the oxygen rich Mira IK Tau 
showed that these carbon bearing molecules can efficiently 
form in the inner regions close to the stellar photosphere 
(DCW99). In this case, it was shown that applications 
of stellar pulsation induced shocks in a narrow region of 
the photosphere can give rise to molecular processes in 
the immediate cooling layer and the hydrodynamic cool- 
ing part (= excursion) of the post shock region and can 
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produce certain carbon bearing molecules including HCN 
and C0 2 . 

It was earlier noted for S-type stars (Bujarrabal et 
al. 1994) that their molecular line strengths are almost 
equidistant between those of C- and O- rich objects. This 
property has been argued to be related to the fact that 
the atmospheric C/O ratio is close to the value of 1 and 
intermediate to that of O- and C- rich evolved stars. We 
therefore apply DCW99's chemical model, which success- 
fully reproduces the observed HCN abundances in O-rich 
objects, to the intermediate case of S stars. In tandem we 
have looked for observational evidence that the HCN is 
forming in the inner envelope the S star \ Cyg. Although 
HCN has previously been detected in \ Cyg, the lines ob- 
served are dominated by material at large radii, so here 
we look for high excitation transitions. In carbon stars, 
such transitions are detected in the inner regions of the 
winds (eg. Lucas & Guilloteau 1992). 

In Sect. ^ we describe the observations and in Sect. |^ 
demonstrate that the HCN detected must lie in the in- 
ner parts of the circumstellar envelope, a result which is 
strongly supported by our modelling of the inner wind 
region of x Cyg (Sect. [|). Our conclusions are in Sect. §. 

2. Observations and results 

Observations were made during 1999 using the James 
Clerk Maxwell Telescope (JCMT) on Mauna Kea, Hawaii. 
We observed the vibrationally excited HCN (0, l lc , 0) J = 
3-2 line at 266 GHz, the (0, l ld , 0) and (0, 2°, 1) J = 3-2 
lines at 267 GHz, and HCN (0, 0, 0) J = 8-7 at 709 GHz. 
The pointing centre was the Hipparcos position of x Cyg, 
19 h 50 m 33. s 92 +32° 54' 50."6 (Perryman et al. 1997). Ob- 
serving parameters are given in Table [l]. 

We detected HCN J = 8-7 and the (0, l lc , 0) J = 3-2 
line at 266 GHz. The spectra are shown in Fig. |l|. The 
HCN (0, l ld ,0) and (0,2°,0) lines at 267 GHz were not 
detected: though this (.'2 = 1 line should be the same 
strength as the one detected at 266 GHz, this spectrum 
had a higher noise level (the frequencies are marked in 
Fig. 0). The HCN and H 13 CN J = 3-2 ground state lines 
at 266 and 259 GHz were also detected. 
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Table 1. Observing parameters: transitions and frequen- 
cies; receivers; beam diameter (full width half maximum); 
zenith opacity at 225 GHz; main beam efficiency and sys- 
tem temperature; RMS noise levels on 1 km s" 1 channels. 



line 


HCN 3-2 


H 13 CN 3-2 


HCN 8-7 




(0,0,0)&(0,l lc ,0) 


(0,0,0) 


(0,0,0) 


v/GRz 


265.8862 & 


259.0118 


708.8772 




265.8527 






receiver 


A3 


A3 


W 


FWHM 


18" 


18" 


8" 


T225 


0.04 


0.1 


0.04 


Vmb 


0.69 


0.69 


0.30 


T 

sys 


270 


400 


3500-5000 


Tmb (rms) 


0.014 


0.030 


0.28 



Fig. 1. HCN spectra towards x Cyg: (top) HCN 8-7; (cen- 
tre) HCN (0, 0, 0) and (0, l lc , 0) J - 3-2 (also offset and 
scaled by 1/100); (bottom) H 13 CN 3-2 with the positions 
of the (0, l ld , 0) and (0, 2°, 0) non-detections indicated. 
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3. Excitation analysis 

The excitation requirements of the J = 8-7 and 
(0, l lc ,0)J = 3-2 transitions are high. We consider 
whether collisional or radiative excitation mechanisms can 
produce the observed line intensities. 



The column densities in each state can be determined 
from the observed integrated line intensities, assuming 
they are optically thin. For HCN this conversion is: 

Vcm- 2 y V GHz / ( J + !) V Kkms" 1 / 

where J is the rotational quantum number of the transi- 
tion J + 1 — *- J, and we assume excitation temperatures 
much higher than the 2.7 K background. The resulting 
column densities are given in Table ^. 

We use the H 13 CN (0,0,0) J = 3-2 as a surrogate 
for the equivalent H 12 CN line as it is more likely to be 
optically thin. We assume the same excitation in H 13 CN 
as H 12 CN, and a [ 12 C]/[ 13 C] ratio, which is unknown in 
X Cyg. Sopka et al. (1989) found [ 12 C]/[ 13 C]> 5 from 
H 12 CN/H 13 CN 1-0 observations. From the ratio of our 
H 12 CN and H 13 CN transitions, [ 12 C]/[ 13 C] > 12. (This 
ratio also supports our assumption that the H 13 CN is op- 
tically thin.) However, these low limits do not rule out a 
much higher value. We take as an upper limit a typical 
ISM value of 70. 

Assuming purely collisional excitation, 

^^M^ + i)^. (2) 

where N\j, ATl, 5u and c/l are the upper and lower state 
column densities and degeneracies and njj 2 the unknown 
hydrogen column density. The critical densities for both 
transitions are high: for the (0, l lc , 0) J = 3 state we calcu- 
late n crit ^ 4.8 x 10 11 cm" 3 (following Stutzki et al. 1988), 
and for J = 8-7 n crit ~ 1.2 x 10 9 cm~ 3 . Applying Equa- 
tion H together with the column densities from Table |[ 
the minimum density requirements to achieve the observed 
ratios are 1 x 10 9 and 5 x 10 7 cm" 3 for (0, l lc , 0) J = 3-2 
and J = 8-7 transitions respectively. For J = 8-7 we as- 
sume the population in J — 7 is the same or less than that 
measured in J = 3. 

The envelope gas density falls off very rapidly with 
radius (Bertschinger & Chevalier 1985; Cherchneff et 
al. 1992). Assuming the stellar parameters for x Cyg given 
in Table 0, the total hydrogen density at IR± is 3.6 x 
10 15 cm -3 and the density falls below the density require- 
ment for (0, l lc , 0) J = 3-2 within 3.5i?* and for J = 8-7 
within 5-R*. Alternatively, assuming a density profile due 
to steady mass loss of 1.8 x 10~ 7 Mq yr" 1 at 8.9 km s" 1 
(BL94), then the density at li?* is 7.5 x 10 8 cm~ 3 and falls 
as R~ 2 . This constrains the HCN 8-7 emission to within 
~ 4i?*, assuming it is collisionally excited. 

An alternative means of exciting these HCN transi- 
tions is by radiative pumping by the IR radiation from 
the star. The (0, 0, 0) and (0, 1, 0) states are connected by 
14/xm radiation with an Einstein A coefficient of 3.7 s _1 
(Ziurys 1986). If collisions can be neglected, and assuming 
geometrical dilution of the stellar radiation field, the ratio 
of column densities in the ground and %)% = 1 vibrational 
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Table 2. Integrated intensities and upper state column densities for detected transitions. 



transition 


HCN 3-2 v 2 = 


H 13 CN 3-2 v 2 = 


HCN 3-2 v 2 = 1 


HCN 8-7 


frequency /GHz 


265.8862 


259.0118 


265.8527 


708.8772 


/ T MB dv/K km s _1 
iVy/cm" 2 


20.31 


1.38 


0.15 


0.30 


4.8 x 10 12 


3.4 x 10 11 


3.6 x 10 10 


1.1 x 10 12 



states depends on the distance from the star according to 
the formula: 

= (eWfc3 ; _ 1){2R/R ^2 + 1 (3) 

with R radius and stellar temperature. Assuming a 
stellar temperature of 2200 K, this constrains the radius 
at which the (0, l lc , 0) line is emitted to R < 17i? + . 

The (0,0,0) J = 8 state cannot be directly pumped 
by radiation (as A J = 0, ±1 for HCN) but it could be 
populated through successive excitations to the (0,1,0) 
levels in increasing J states. The radius requirements for 
this to take place are stronger than for the (0, l lc , 0) state. 

For either collisional or radiative excitation, the ra- 
dius requirements are tightened further if either (a) 
[ 12 C]/[ 13 C]< 70 or (b) H 13 CN J = 3-2 emission origi- 
nates from a larger region than HCN (0, l lc , 0) J — 3-2 or 
HCN J = 8-7. 

To conclude, the excitation requirements are such that 
both the HCN (0, l lc , 0) J = 3-2 and HCN J = 8-7 tran- 
sitions must originate from gas within < 20 stellar radii, 
in the inner envelope of the star, whether collisionally or 
radiatively excited. 

4. Model of the inner wind of x Cyg 

We have used our chemical model, which has succesfully 
produced the observed HCN abundance in the case of an 
O-rich object (see DCW99 for details), with a value of 
C/0 ~ 1 for modelling the inner wind region of x Cyg. 
The stellar parameters considered in this study are listed 
in Table [|. The radius is obtained from the observational 
value of Tuthill et al. (1999) and the pulsation period is 
from Bedding & Zijlstra (1998). Using the standard pul- 
sation equation for Miras of Fox and Wood (1982) and as- 
suming that the star is pulsating in its fundamental mode 
(Q =0.09), we derive a stellar mass of 1.15 Mq, which is 
a typical value for the stars of this class. The temperature 
of 2200 K was obtained from Haniff et al. (1995) which 
gave a luminosity which is close to the canonical value for 
typical Miras. 

We have considered here the inner wind above the pho- 
tospheric region which experiences passage of strong, pe- 
riodic shocks generated by stellar pulsation. The model 
deals with the chemistry of the immediate region (thermal 
cooling region) and the hydrodynamical cooling region of 
post shock described by Bertschinger & Chevalier (1985), 
Fox & Wood (1985) and Willacy & Cherchneff (1998). 



Table 3. \ Cyg - stellar parameters: 7 and a are defined 
as in Cherchneff et al. (1992) and X (Y) = X x 10 Y . 



D 106 pc M l-8(-7) M yr" 1 

2200 K L* 1.8(3) L Q 

fU 290 R Q M* 1.15 M 

P 408 days n(r shock ) 3.6 (15) cm" 3 

r shock 1.0 R* a 0.6 



7 = fshock/fesc 0-89 C/O ratio 0.95 



Table 5. Calculated fractional abundances (relative to 
the total gas number density) versus shock strength and 
radius. 



Species 


T.E. 


32 km/s 


26.1 km/s 


22.6 km/s 




1. R* 


1. R* 


1.5 R* 


2. R* 


H 


1.0(-01) 


4.3(-03) 


4.4(-03) 


4.2( 


-01) 


H 2 


7.3(-01) 


8.2(-01) 


8.1(-01) 


4.4( 


-01) 


S 


1.5(-05) 


2.0(-06) 


2.9(-06) 


1.6( 


-05) 


C2H2 


1.5(-15) 


3.6(-ll) 


4.6(-12) 


4.9( 


-15) 


cs 


3.4(-09) 


1.5(-05) 


2.0(-05) 


3.0( 


-06) 


HS 


4.8(-06) 


4.3(-06) 


l.l(-06) 


7.0( 


-09) 


H 2 S 


4.8(-07) 


4.8(-07) 


4.8(-08) 


7.2( 


-13) 


NS 


4.6(-ll) 


7.3(-12) 


2.9(-ll) 


9.4( 


-13) 


N 


3.4(-09) 


3.0(-13) 


7.6(-17) 


1.1( 


-17) 


N 2 


9.0(-05) 


5.8(-06) 


5.3(-05) 


7.1( 


-05) 


HCN 1 " 


1.6(-09) 


1.8(-04) 


8.5(-05) 


4.8 


(-06) 


CN 


1.4(-11) 


4.6(-09) 


6.2(-10) 


1.9( 


-09) 


O 


2.6(-08) 


8.4(-ll) 


l.l(-ll) 


8.7( 


-10) 


OH 


2.0(-07) 


3.0(-08) 


3.7(-09) 


1.1( 


-08) 


H 2 


7.4(-06) 


2.0(-04) 


l.l(-04) 


5.6( 


-06) 


CO 


9.9(-04) 


8.5(-04) 


8.6(-03) 


6.6( 


-03) 


COa 


2.0(-09) 


6.4(-08) 


8.3(-07) 


4.7( 


-06) 


2 


1.2(-13) 


6.8(-15) 


5.5(-16) 


5.3( 


-15) 


SiO 


4.4(-05) 


4.7(-05) 


5.0(-05) 


3.9( 


-05) 



The observation value for HCN ( averaged over the entire 
envelope) is 2.5(-6) derived from millimeter observations 
(Bujarrabal et al. 1994) 



4-.1. The chemistry 

We have considered 760 reactions involving 68 chemical 
species and have all possible chemical routes possible in 
a dense gas. The reaction rates considered are the same 
as in DCW99 (and the references therein). We have as- 
sumed thermal equlibrium for the photosphere and de- 
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Table 4. Pre-shock, shock front and excursion (= post-shock) gas temperature and number density as a function of 
position in the envelope and shock strengths. M is the Mach number associated with each shock speed. 



Position 


Shock Vel. 


M 


Pre-shock 


Shock Front 


Start of 


excursion 


(**) 


(km s" 1 ) 




To (K) n (cm" 3 ) 


T (K) n (cm" 3 ) 


T (K) 


n (cm -3 ) 


1.0 


32.0 


10.3 


2200 3.62 (15) 


47458 2.07 (16) 


6528 


1.26 (17) 


1.5 


26.1 


9.51 


1724 1.28 (13) 


31959 7.28 (13) 


4867 


3.98 (14) 


2.0 


22.6 


8.98 


1452 3.93 (11) 


24128 2.22 (12) 


3966 


1.13 (13) 



rived molecular abundances for the temperature, C/O ra- 
tio and gas density given in Table 3. We then apply a 
shock to the photosphere (the shock velocity considered is 
Vshock— 32 Km s _1 ) and study the chemistry in the im- 
mediate and excursion region. The resulting abundances 
form the input for the shock at the next distance. The 
different parameters describing the shock structure at dif- 
ferent radial distances from the star is given in Table 4. 

4-2. Model results and discussion 

The molecular abundances of few selected species relative 
to the total gas number density are given in Table | for 
various shock strengths. We have considered shock chem- 
istry at gas layers very close to the star. The abundances 
at 2R± may not be the exact values in the outflow because 
of uncertainties in the dust formation radius which halts 
the shock chemistry, but are indicative. One can see that 
in addition to molecular hydrogen certain other molecules 
like CO, H2O and N2 are dominant, confirming that they 
are the parent molecules. But, to come to the point of 
our interest, the shock chemistry is responsible for forma- 
tion of species like HCN along with CO2 and CS close 
to the star. The prediction of HCN at the distance and 
environment described is in excellent agreement with the 
observational result of Section ||. Moreover, the theoreti- 
cal value obtained for the HCN abundance is in agreement 
with that derived from millimeter line observations in the 
outer wind. 

The chemical processes responsible for the formation 
of HCN involves reactions with cyanogen - an observation 
which was earlier found to be the case in O-rich IK Tau 
(DCW99). HCN is formed by the reaction 

CN + H 2 -> HCN + H (4) 

CN acts as an intermediary in the formation of HCN and 
CS and is destroyed quickly by atomic hydrogen. The re- 
action is sensitive to temperature and is very fast in the 
gas excursions. On the basis of the chemistry we can claim 
that HCN is a direct result of the shock chemistry in the 
inner wind region that we have considered here and travels 
as a parent species through the envelope unaltered until it 
reaches the photo-dissociation regions in the outer wind. 
The choice of the ratio C/O=0.95 is arbitrary, but with 
different values of C/O equal or slightly greater than 1.0, 
though the T.E. abundance values differ by a factor of 100, 



yet the resulting HCN abundance at 2R* remains within 
5 percent of the quoted number in Table 5. 

5. Conclusions 

The detection of HCN in high excitation states clearly 
demonstrates that it is formed in the inner part of the 
envelope of x Cyg, within < 20-R*. This is strongly sup- 
ported by detailed chemical modelling of the atmosphere. 
The earlier claim of HCN along with other carbon bear- 
ing molecules being formed in the innermost region of the 
envelope of O -rich star like IK Tau (DCW99) seems to 
be true in the case of S type star x Cyg as well. This is 
the first observational verification of the theoretical claim 
that shock chemistry plays an active role in defining the 
chemical composition of the inner wind regions of Miras. 
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